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The Gluon and Ghost propagators in Landau gauge QCD are investigated using the Schwinger- 
Dyson equation approach. Working in Euclidean spacetime, we solve for these propagators using 
a selection of vertex inputs, initially for the ghost equation alone and then for both propagators 
simultaneously. The results are shown to be highly sensitive to the choices of vertices. We favor 
the infrared finite ghost solution from studying the ghost equation alone where we argue for a 
specific unique solution. In order to solve this simultaneously with the gluon using a dressed-one- 
loop truncation, we find that a non-trivial full ghost-gluon vertex is required in the vanishing gluon 
momentum limit. The self-consistent solutions we obtain correspond to having a mass-like term in 
the gluon propagator dressing, in agreement with similar studies supporting the long held proposal 
of Cornwall. 
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I. INTRODUCTION 

The Gluon Schwinger-Dyson equation and its gauge 
fixing counterpart, the Ghost Schwinger-Dyson equation, 
are principal tools in investigating the theoretical prop- 
erties of strongly coupled QCD. The transition from the 
weak coupling regime at large momenta, where pertur- 
bation theory applies, down to strong coupling at small 
momenta gives rise to the key non-perturbative phe- 
nomena of confinement and dynamical chiral symmetry 
breaking. The latter endows quarks, and perhaps glu- 
ons, with mass. These are not physical pole masses, but 
rather effective Euclidean masses dynamically generated 
by non-perturbative interactions [1, 2]. The main reason 
for studying these propagators is that they are a nec- 
essary input for the Bethe-Salpeter (BSE) and Faddeev 
equations used for calculating physical quantities, such 
as hadron masses and form- factors [3, 4]. 

The infinite tower of Schwinger-Dyson equations 
(SDEs), and its necessary truncations in QCD, can be 
a troublesome beast. This paper is intended to stimulate 
discussion about the assumptions commonly made and 
the solutions to which these lead. We believe that im- 
portant issues remain unaddressed and these are outlined 
below. 

The solutions provided by the SDEs may be compared 
with the results of complementary techniques, such as 
Lattice QCD. A benefit of this is that the theoretical 
strengths and weaknesses of each method fall in different 
areas. Both may be formulated in precisely the same way, 
using the same gauge and freely varying the quark con- 
tent. While lattice calculations are inevitably restricted 
to larger quark mass, an advantage of the SDE method 
is that there are no major computational obstacles to ex- 
tending these to the physical light quarks with masses of 
O(MeV). 

The paper is organised as follows, in Sec. 2, we intro- 



duce the Schwinger-Dyson equations for the gluon and 
ghost system and proceed to give details regarding trun- 
cations. In Sec. 3 we solve the ghost equation for simple 
truncations and show that solutions may be robustly ob- 
tained for a wide range of vertex inputs. In Sec. 4 we out- 
line the gluon equation and some of the issues we deem 
important that deserve immediate attention. We present 
simultaneous solutions for the gluon and ghost propaga- 
tors for a range of assumptions to indicate the level of 
precision currently available. In Sec. 5 we compare our 
results to those obtained using Lattice QCD simulations. 
Section 6 briefly gives our conclusions and outlook. 



II. SCHWINGER-DYSON EQUATIONS FOR 
GLUONS AND GHOSTS 

The full equation for the gluon propagator, in the ab- 
sence of quarks, is represented in Fig. 1. Without ap- 
proximations, this is given by, 

Or 1 = Df^\ P y l + ng M „(p) + 

+ K ,u (p) + K (p) + K ^(p), (i) 

where latin indices denote colour and greek letters de- 
note the Lorentz suffixes. The loop integrations II(p) 
are labeled by the number of particles in the loop with 
c corresponding to ghosts and g corresponding to glu- 
ons. Eq. (1) is represented graphically in Fig. 1. The 
momenta flowing through the internal propagators are 
arbitrary. However we use a symmetric routing for the 
one-loop diagrams with l± — £±p/2, where i is the loop 
integration momentum and p is the external momentum 
that is divided equally through each propagator. When 
working with a large finite cutoff as is often done numer- 
ically, then a symmetric arrangement is often essential to 
preserve translation invariance. 
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FIG. 1. The Schwinger-Dyson Equation for the Gluon in the 
absence of quarks. 

We work in covariant gauges. Then neglecting quarks, 
the gluon equation only couples to ghosts, and these have 
their own Schwinger-Dyson equation, which is depicted 
in Fig. 2 and given by, 

oW^^^ + n^), (2) 

where H gc contains just one simple loop integration. 

Solving these equations depends non-trivially on the 
vertices, which in turn depend on higher n-point Green's 
functions and therein lies the key issue. Truncation of 
this sequence of dependent Green's functions is necessary 
to be able to make predictions. However truncating has 
the potential to introduce errors and violate the physical 
properties of the theory. 

A. Sketching the general features of the equations 

The interplay of the various terms in the equations 
and identifying which terms are important will be our 
primary concern. For the most part we will adopt a trun- 
cation in which only one loop dressing of the gluon prop- 
agator will be considered. This is a favorite truncation 
scheme [5-11]. There are then four quantities of inter- 
est seen in Fig. 1: the two propagator dressings that the 
equations should determine, and the two unknown ver- 
tices that are required inputs. There are important inter- 
relations between these terms constrained by the struc- 
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FIG. 2. The Schwinger-Dyson Equation for the Ghost. 



ture of the SDEs and the Ward-Slavnov- Taylor identities 
(WSTIs). 

The propagators are remarkably simple. The ghost 
and gluon dressings are just functions of p 2 with the gluon 
dressing multiplied by a tensor structure transverse to its 
momentum p. In the Landau gauge these are, 

K(P 2 ) = ^ bG f 1 (^- P ^) (3) 

V ab (p 2 ) = S abGh{ ^ (4) 

where V^ v is the full gluon propagator, V is the ghost 
propagator. The a, b indices relate to the color carried 
by each propagator. The functions Q£ and Qh are the 
respective dressing functions, for which we will be solving 
later. These functions contain all of the non-perturbative 
physics of these two Green's functions. 

The two propagator dressing functions are determined 
from the various quantities that make up their equations, 
which are represented diagrammatically in Fig. 1 and 
Fig. 2. The triple-gluon vertex only appears in the gluon 
equation. However the two propagators and the ghost- 
gluon vertex appear in both equations. While the dress- 
ing functions are not expected to change sign over the 
whole Euclidean momentum region, the different contri- 
butions from the loop diagrams can have different signs, 
and indeed they do. The ghost-loop has an additional 
( — 1) owing to the antisymmetry of the ghost field. 

The gluon propagator dressing determined in the mo- 
mentum subtraction scheme is given by, 

^l-'^Vl^ + n^V^ + n^pV 2 ) (5) 

where p, 2 is the subtractive-renormalization point, 
^2c(p 2 ,p 2 ) is the ghost-loop contribution subtracted 
from itself at fi 2 , similarly for H2 g (p 2 , /x 2 ) which is the 
gluon loop contribution. The precise details may be 
found in the appendices. 

In the definitions we use, the gluon propagator dress- 
ing function Ql(p 2 ) is positive everywhere. However, the 
dressed-gluon-loop diagram gives a negative contribution 
in the momentum subtraction scheme, below the sub- 
traction point. This is important because it tells us that 
the two other terms on the right of Eq. (5), Q£{^ 2 ) and 
^2c(p 2 , M 2 ) must add to give an overall positive contribu- 
tion of greater magnitude than that of Yl2 g (p 2 , p 2 )- Now 
in the perturbative region Q£(pi 2 ) — 1 is more than large 
enough to ensure that Ql{p 2 ) remains positive, but as 
we evolve down to infrared (IR) momenta, then the two 
one-loop contributions are both large. Quite generically 
for many sensible vertices, the ghost-loop term must then 
give a large positive contribution in order to be able to 
solve the equation. A dichotomy of explanations exist: 
the large contribution either arises from a singular ghost 
propagator, or a non-trivial ghost-gluon vertex. Despite 
this distinction, little change is observed in the gluon 
dressing function since the overall contribution from the 
ghost loop remains similar. 
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The full ghost-gluon vertex has two structures that 
may be written as follows, 

r ii (k,p,q) = igf ahc [q lt a(k,p,q) + k^(k,p,q)} , (6) 

where the momenta arc defined in Fig. 3 and a is the non- 
perturbative function dressing the outgoing ghost mo- 
mentum and /3 is the non-perturbative function dressing 
the outgoing gluon term. This can be arranged in a num- 
ber of ways since the momenta are related by k = p — q. 
In the ghost equation the term j3 does not contribute in 
Landau gauge, however in a properly projected gluon the 
j3 term can be important. It appears in the ghost-loop, 
as we emphasise later. 

The ghost-gluon scattering kernel is an important 
quantity since it appears in dressing both vertices. It is 
directly related to the ghost-gluon vertex, 

T^k,p,q)=igf abc q v t^(k,p,q) . (7) 

Hence the bare version is T M „ = g^ v in order to yield 
the tree-level vertex. The same function also appears in 
the WSTI for the triple-gluon vertex, which is a feature 
of the interconnections between terms in a gauge theory. 
The precise relation will probably be necessary to remove 
all gauge-dependence from the results. We do not give 
further details here as these are fully covered in [6, 12] 
and references therein. 



B. Gauge fixing 

We work in Landau gauge since it happens to be the 
most theoretically appealing for several reasons. An obvi- 
ous and oft-suggested alternative is to use a gauge defined 
by an axial vector n M , where ghosts do not appear. How- 
ever these have other drawbacks, not the least of which 
is the breaking of Lorentz symmetry, and in particular 
in numerical studies where l/(n.p) terms are a complica- 
tion. Additionally, direct comparison to Lattice QCD is 
less straightforward in these gauges. The radiative cor- 
rections to the Gluon propagator are always transverse in 
the external momentum due to the Slavnov- Taylor iden- 
tity, and in the Landau gauge the Gluon itself is also 
transverse, which is a simplifying feature. Furthermore, 
there are theoretical statements regarding confinement, 
gauge fixing and the vanishing momentum limit of QCD. 
These are usually formulated in Landau gauge and some 
details are given below. 

Ultimately we hope that these methods will be useful 
in making physical predictions where the gauge depen- 
dence must drop out as in perturbation theory. However, 
the gluon propagator dressing function is a gauge depen- 
dent object and one expects it will combine with all of 
the other gauge dependent functions in making gauge- 
independent physical observables. 

A deeper analysis of the gauge fixing procedure is re- 
quired on theoretical grounds, since the Faddeev-Popov 



procedure may not be sufficient non-perturbatively. Ear- 
lier studies by Gribov, Kugo, Ojima, and Zwanziger have 
received refinement in recent years in light of the serious 
debate over the behavior of the propagators in the van- 
ishing momentum limit [13-16]. Formulated in Landau 
gauge, the Gribov-Zwanziger gauge fixing procedure con- 
siders the technical issues related to gauge field Gribov 
copies and restriction of the gauge field space to select 
one representative configuration. The broad conclusion 
at present is that the dressed gluon propagator carrying 
momentum p is suppressed as p vanishes, with respect 
to its bare counterpart. Moreover, many studies find 
a propagator that goes to a constant, or equivalently a 
propagator dressing function that is proportional to p 2 
in the small p 2 limit. 

Coupled with the earlier conclusions of Kugo and 
Ojima [17] regarding a singular ghost dressing being a 
signal of confinement, SDE practioners were led to search 
for a particular class of solution where the leading powers 
in the IR of the ghost and gluon propagators are directly 
linked. More details of this class of solution will be out- 
lined below, although the key feature is the singular ghost 
dressing function. Recently, new studies combining these 
two earlier studies, and also additional refinements using 
Stochastic quantization techniques, have found that a fi- 
nite ghost dressing solution is certainly allowed and may 
even be preferred [18-22]. 



C. Existing Truncation Schemes and Solutions 

We now turn to the act of truncation and briefly sum- 
marise some of the ideas previously considered for these 
equations. Most truncation schemes for solving the gluon 
equation are motivated by the practicalities of finding 
solutions. For example, in all but one of the references 
in the literature [23] the dressed-two-loop diagrams of 
Fig. 1 are entirely neglected. This is not due to some 
convenient power counting scheme, but rather because 
their inclusion is both numerically challenging and com- 
putationally intensive. Clearly in the perturbative ultra- 
violet (UV) regime, these terms are subleading. However 
in the intermediate and vanishing momentum limit, their 
contribution is essentially unknown. Naively, one might 
expect that if the gluon is heavily suppressed in the small 
momentum limit then these terms may be unimportant. 
In this case of an IR suppressed gluon, the most impor- 
tant region for these terms would be the mid-momentum 
region (i.e. p ~ 1 GeV) as the leading perturbative log- 
arithm reaches a maximum and turns over in the non- 
perturbative region. For a gluon propagator that is con- 
stant in the IR (or equivalently a gluon propagator dress- 
ing function that goes as p 2 ), then the situation is less 
clear. The answer depends on the vertex ansatze adopted 
and the full solution of the gluon dressing function. The 
mid-momentum region is also the most physically rele- 
vant so eventually these terms must be properly included. 

A recent one-loop perturbative study, in which an effec- 
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tive gluon mass term is artificially introduced, has found 
good agreement [24] with results from Lattice QCD. This 
differs from a typical SDE study since the mass is inserted 
by hand and iterations are not performed. Otherwise the 
equations are very closely related. In order to perform 
iterations required for self-consistent solutions (see Ap- 
pendix A), some term in the gluon equation is required to 
dynamically generate the mass term. In the SDEs, this 
type of contribution could come from any of the loops in 
the gluon equation including the two loop graphs. 

Nevertheless, we proceed by neglecting these two-loop 
dressing terms and turn to the remaining two unknown 
vertices that appear in the ghost and gluon equations. 
The ideal situation would be that the propagators de- 
pend only weakly on the ansatze for these vertices, and 
it is the propagator itself and the structure of the SDE 
that determines the result. In studies in QED the situa- 
tion is indeed close to this ideal, since the Ward-Green- 
Takahashi identity determines the key vertex dressings in 
terms of inverse propagators [25-30]. QCD is inevitably 
more complicated. The Ward-Slavnov- Taylor Identities 
that relate vertices and propagators do not admit such 
simple solutions as in QED, and for the full vertex there 
is a complex interplay between the gluon and ghost sec- 
tors. 

^ t k 
O 

o 
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FIG. 3. The ghost-gluon vertex indicating the momentum 
definition we adopt, the outgoing ghost momentum is q and 
the gluon momentum is k. 

In QCD the starting point has often been bare ver- 
tices, or vertices dressed by simple ratios of dressing func- 
tions [3, 5, 8, 31-34]. In particular, one option is to use 
the bare ghost-gluon vertex in place of the full vertex, 

rf c (k lP ,q) -> Tf abc {k,p,q) = igf ahc % , (8) 

which at first sight seems like an over-simplification, since 
this cannot be correct even at the perturbative one-loop 
order [35, 36], let alone non-perturbatively. However, it 
is expected that the ghost-gluon vertex should be rela- 
tively simple in the Landau gauge, by appealing to Tay- 
lor's non-renormalization theorem [37] which is valid in 
carefully chosen renormalization schemes [38, 39]. The 
full vertex is also expected to reduce to its bare form for 
vanishing incoming ghost momenta [7, 37, 38], 

^mT^p^) = T^ abc (k,p,q), (9) 

the momentum definition used throughout being indi- 
cated in Fig. 3. This simple vertex appears to work for 



the infinite IR ghost dressing function [8, 23]. However in 
the finite solutions a problem arises. In practical terms, 
the large contribution from the infinite ghost must be re- 
placed by a large contribution from the vertices in order 
for the gluon dressing function to remain positive, and 
for self-consistent solutions to be obtained. There are 
two forms given in the literature that accomplish this. 
Firstly, a vertex that is transverse to the gluon momen- 
tum in the IR is suggested as a possible solution [7, 11]. 
This is motivated by solutions in which the ghost is dom- 
inant in the IR region. Then the ghost loop controls cor- 
rections to the gluon propagator, see Fig. 1, which must 
then be transverse on its own in this regime. This can 
be achieved by making the ghost-gluon vertex transverse 
to the gluon momentum in the IR. A form that achieves 
this is [7, 11], 

r« - igf abc - k^T m (k lP ,q^j . (10) 

where J-jr is a smoothed step function defined in 
Eq. (B15) that switches on this behavior in the IR, but 
ensures it does not affect the UV. J-'m, = 1 in the IR 
and must vanish in the UV in order to reproduce the 
perturbative results. To achieve this, the function 
contains a free parameter that controls the momentum 
of the switchover point. This has to be fixed and some 
sensitivity to its value will be present in the solutions. 
This extra term has no effect in the ghost equation, where 
the transverse gluon is contracted with this vertex, and 
so it automatically drops out. It can however appear in 
the gluon equation as we demonstrate below. The super- 
script label (1) in Eq. (10) will be used to refer to this 
vertex later. 

The second form considered in the literature involves 
inserting massless poles into the vertex that give an IR 
enhancement [40-42]. This is in an alternative trunca- 
tion where the Feynman diagrams are collected into in- 
dividually transverse groups using pinch-technique and 
background field method rearrangements [10, 43]. There 
is not a direct comparison of such a vertex insertion in 
this formulation with the standard SDEs we use. We can 
however make similar considerations and investigate the 
effects of massless poles in the vertex. 

Several different choices are made for the triple-gluon 
vertex in the literature, which is perhaps a little supris- 
ing since sensible WSTI solutions are available. How- 
ever these depend on an unknown contribution from the 
ghost-gluon vertex. When solved self-consistently, the 
SDE solutions should be matched on to the resummed 
leading-logarithm from perturbation theory. However us- 
ing symmetric triple-gluon vertices, it has been found 
that the anomalous dimension of this logarithm is not 
exactly reproduced [5]. 

Proposed solutions to this issue include solving an 
equation relating the vertex dressing and the anomalous 
dimension in the perturbative region in order to choose a 
vertex dressing that gives the correct outcome. Typically 
these vertices sacrifice Bose-symmetry between the inter- 
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nally contracted and external legs of the vertex. When 
a dressed symmetric vertex is considered, this diagram 
then contributes in the IR region of the gluon equation. 
This in turn requires an interplay between the ghost and 
gluon loops in the IR for the gluon equation and ar- 
guments about an individually transverse ghost-diagram 
may not necessarily be valid [7, 11]. 
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FIG. 4. A sketch of typical dressing functions. The momen- 
tum scale is arbitrary, but can be thought of as ~ 1 GeV. 

There are two broad classes of solution in existence 
and these have been strongly debated [10, 11, 44]. The 
gluon, shown in Fig. 4, is qualitatively similar in most 
recent studies. However the ghost in the IR can vary 
wildly: from a finite value to infinity. We consider the 
ghost equation in detail in the next section. There are 
also solutions from Lattice QCD that are produced in 
the absence of quarks. At present these point to a fi- 
nite ghost propagator dressing function [45]. We give 
an example of some solutions in Fig. 4 including a typi- 
cal gluon dressing with the logarithm from perturbation 
theory at large p 2 and a suppression at small p 2 . The 
ghost is less certain and a wide range of solutions may 
be found. The most important region is around 1 GeV, 
where QCD becomes confining. How well the functions 
can be determined in this region will be a key concern. 
The vanishing momentum region is physically less rele- 
vant, although the interplay between the ghost and gluon 
there does determine whether or not self-consistent solu- 
tions can be obtained. Importantly, the method we use 
to find such self-consistent solutions is outlined in Ap- 
pendix A. 

III. THE GHOST EQUATION 

Non-linear integral equations may admit multiple solu- 
tions and it is apparent that the ghost propagator equa- 
tion falls into this class. Early attempts to include ghosts 
in the solution of the gluon equation used the well known 
infrared (IR) power law assumption for the vanishing p 2 



limit of the propagators [46, 47], 

Um ge(p 2 ) = a g ( P 2 r*, (11) 

p 2 — >0 

lim Gh(p 2 ) = a c {p 2 T% (12) 
p 2 — >o 

the relation that is found is K g — 2k and k c = —k. 
This behavior, where the power k is linked, is derived 
from matching powers on both sides of the ghost equa- 
tion using a bare ghost-gluon vertex. In solving the ghost 
equation the coefficients a g and a c are left unconstrained. 
However the gluon equation fixes a relation between the 
two. Typically, vanishing IR gluon propagators are then 
considered and hence a singular ghost arises for this spe- 
cific solution. 

More recent studies [11, 38, 48-50] have found a second 
solution where the IR powers are not linked; a vanishing 
gluon admits a finite ghost solution also. There the pow- 
ers are typically n g — 1 and n c — 0. The solution selected 
by physical QCD is still an open question however and at 
present there are ideas about the singular ghost solution 
being a critical endpoint of a family of solutions including 
the finite ones [51], and also suggestions of an additional 
gauge parameter [44, 52]. 

It is possible that solutions in the deep IR are not 
solutions for the whole momentum region and this will 
become important when we consider the gluon equation 
itself. However, both of these classes of solutions can be 
found using numerical methods in Euclidean space using 
standard techniques. There is one key distinguishing fea- 
ture between the two however. The infinite IR ghost so- 
lution requires the ghost to be subtracted (renormalised) 
at zero momentum. Its value for physically relevant mo- 
menta is then specified. How the ghost evolves is fixed. 
The finite solution may be found by subtraction at any 
value of momentum. 

We presently do not state a strong preference for ei- 
ther solution. However we note that only the finite solu- 
tion is allowed by a perturbative renormalization scheme 
without fine-tuning. The infinite solution never arises 
without such tuning, one must search for it by specifying 
the infinite value at zero momentum and sacrificing such 
freedom in the UV. We will elaborate on this below. 



A. Solutions of the ghost equation alone 

Using a fixed gluon input we may solve the Ghost SDE 
alone and investigate its sensitivity to a range of input 
vertices and gluons. This is useful because this type of 
SDE, or gap equation, is very simple to solve, and it 
teaches us what to expect when solving the more com- 
plicated, coupled equations self-consistently. 

In doing this we utilise the bare vertex of Eq. (8) , and 
the following model gluon [53], 

SV) = w (13) 

- 2 +P 2 (l + ^^Log(^)) 22 
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FIG. 5. (Color Online) Examples of ghost solutions on a log- 
log plot subtracting at zero momentum. Only the specified 
subtraction value Gh(0) is varied between the solutions which 
can be read off. The Gh(p 2 = 0) value is fixed and the ghost 
equation solved with the depicted gluon until the ghost inputs 
and outputs are self-consistent. The gluon dressing Q£(p 2 ), is 
the dashed curve and the solid curves are the different ghost 
dressings Qh(p 2 ). The units of p 2 are arbitrary since we have 
not fixed the coupling to the physical value, but may be con- 
sidered to be (I GeV 2 ). 

which reproduces the correct UV perturbatively re- 
summed logarithm, see Eq. (B16), and provides an IR 
mass term necessary to produce a gluon propagator 
dressing function that behaves as p 2 in that limit. We 
utilize the parameters p? = 10 4 , m 2 = 0.1, g 2 (n 2 )/4n — 
0.12 and p = 1 which give the gluon shown in Fig. 5. Di- 
mensionful quantities can be thought of as having units 
close to GeV. However we do not match to a phys- 
ical scale at this stage. The physical scale is deter- 
mined by the value of the coupling at the renormaliza- 
tion point. This of course only has true physical mean- 
ing when quarks are included. This gluon, Eq. (13), 
is qualitatively similar to that found in recent Lattice 
QCD studies [45, 54] and in other Schwinger-Dyson stud- 
ies [11, 12, 49, 53, 55]. 

When solving the ghost equation and subtracting at 
zero momentum, we find its value at the subtraction 
point does not change its UV values between many of 
the solutions. Thus, for ghosts with IR values in the 
range 2 < Qh(0) < oo we find their UV differences to be 
negligble. Note that the precise value is dependent upon 
the coupling and the gluon equation. For a fixed gluon, a 
larger coupling gives a larger critical value of Qh(0) above 
which all the dressings are practically indistinguishable 
in the perturbative region. This can be seen in Fig. 5 
where the largest three ghost dressings, although differ- 
ent in the IR, are identical in the UV. Below this critical 
value, the ghost equation admits solutions that differ in 
both the UV and the IR. They still exhibit the same per- 
turbative logarithm although it is hardly visible on this 
scale since its coefficient is reduced when Qh(0) is set to 
be small. Plotted on a linear vertical axis, the effect is 



clearer as can be seen in Fig. 6. 

This is interesting and a point that has not often been 
stressed. If we subtract at some perturbative value and 
impose the perturbative condition Qh(p, 2 ) — 1, and do 
the same for the gluon then none of the curves in Fig. 5 
are reproduced. The solution that connects to the stan- 
dard perturbative solution with the standard momentum 
subtraction condition is unique and is separate from both 
the singular ghost and the zero-momentum subtraction 
finite solutions (except for the single fine-tuned value sat- 
ined by this unique solution). 

Setting Gh([i 2 ) = 1 is not essential as we can in prin- 
ciple renormalise at any point and the renormalization 
group tells us how these differently renormalised solu- 
tions are connected to each other. However, it is not 
possible to run down from the solution that we have ob- 
tained subtracting at a perturbative point and enforcing 
Gh(fi 2 ) = 1 to any general solution subtracted at zero 
momentum. Fine-tuned examples exist but these essen- 
tially predetermine the IR value. 

Starting from the perturbative solution and running 
down into the IR, we would then stay on the finite solu- 
tion for the ghost dressing for all momenta. The expec- 
tation is that asymptotically free QCD at large momen- 
tum transfers is accurately described by the perturbative 
solution, so this is the one we favor. In this example, 
a perturbative subtraction point never admits the infi- 
nite IR ghost solution. This effect is depicted in Fig. 6, 
where the solid curve is the unique solution selected by 
renormalising Qh(p 2 ) — 1 at the same /x 2 point as the 
gluon. Hence, given similar conclusions from other stud- 
ies [11, 38, 48-50, 56], we carry forward this solution to 
investigate the gluon equation. 

In contrast, subtracting at zero momentum and speci- 
fying different values for Qh(0) maps out infinitely many 
other solutions even with the same fixed gluon input. 
Evolving these solutions up in momentum yields different 
curves in the perturbative region. We have verified that 
this effect is also present for fixed gluon inputs that van- 
ish more or less rapidly in the IR, including the k ~ 0.6 
solution [8, 9, 11, 57]. 

B. Other ghost-gluon vertices in the ghost equation 

It has been noted [38] that Taylor's theorem, which 
is often used to infer a bare ghost-gluon vertex, actually 
only places a restriction on the sum of the two functions 
dressing the components of the vertex, Eq. (10). The con- 
dition is in the simplest terms, that any corrections van- 
ish when the incoming ghost-momentum vanishes. This 
motivates the modification to Eq. (10), 

If 5 = igf abc [q, - pJ^M m T m {k,p, qfj , (14) 

which is clearly not transverse to the gluon momentum. 
However, the piece additional to the bare term now van- 
ishes when the incoming ghost momentum p vanishes. 
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FIG. 6. (Color Online) The dotted curves correspond to the 
zero momentum subtracted solutions and the colours match 
Fig. 5. The dashed curve is the gluon input, Q£(p 2 ), and 
the solid curve is the physically relevant solution of the ghost 
equation, Qh(p 2 ). The units of p 2 are arbitrary. 

The non-perturbative normalization parameter Air, is 
also introduced into the modeling. Unlike dressing forms 
that add only to the term, this form may affect the 
solutions of the ghost equation, as we demonstrate below 
in Fig. 7. Adopting the same method as above, we fix 
the gluon using the previous model and parameters, and 
solve the ghost until self-consistent with this input and 
its own equation. 
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FIG. 7. (Color Online) A comparison of two solutions nor- 
malised using the perturbative condition. The (red) dotted 
curve is the ghost solution, Gh(p 2 ), with a bare vertex or 
transverse vertex. The (blue) dashed curve is the ghost dress- 
ing obtained using Eq. (14). The (black) solid curve is the 
fixed gluon dressing input, Q£(p 2 ) from Eq. (13). The units 
of p 2 are arbitrary internal units. 

We find that in principle the ghost equation is depen- 
dent upon the choice of vertex and fairly large corrections 
are possible as found for the example shown in Fig. 7. 



Many other extensions are possible and these exist in the 
literature. Clearly a precise form for this vertex would 
be useful. However at present we must consider this part 
of the uncertainty on the result of the calculation. Cor- 
rections of 20% over the bare vertex appear to be al- 
lowed. Nevertheless we have found the ghost equation to 
be solvable for a wide range of parameters and choices of 
vertices. 



IV. COUPLING GLUONS AND GHOSTS 

After the ease of solving the ghost equation, we now 
turn to the gluon, where matters are very different. 

It is the opinion of the present authors that no com- 
pletely satisfactory solution exists for the gluon propa- 
gator equation at present. Solutions exist, but they are 
highly dependent on arbitrary vertex choices, and dis- 
appointingly this is an issue rarely addressed. The pre- 
dictive power of the equations is often lost through the 
introduction of arbitrary parameters that are not derived 
from the fundamental field theory. The simplest possi- 
ble model just modifies the gluon propagator, such that 
it contains an IR mass term that vanishes in the UV. 
Modeling vertices with a number of parameters and a 
range of possible forms is unfortunately no more predic- 
tive. However, it can be useful in guiding where to look 
and in understanding those quantities we may need to 
know precisely and those that are less important. 

Although we intend to subtract and renormalise at a 
perturbative point, which excludes the singular ghost so- 
lution, it is worth saying a few words about one partic- 
ular singular solution, especially in the light of recent 
criticism [57]. In what follows, we note that this is not 
directed at all singular solutions, just this particular one 
and any others that contain the same flaw we describe. 

The problem arises due to the combined effects of trun- 
cation and projection of the dressed gluon propagator in 
the Landau gauge, 

zV(p) = —p- [9^ --pT)> ( 15 ) 

= Atf) 9ltu -B(j?)^. (16) 

The structure of the gauge theory demands that A(p 2 ) = 
Qi{p 2 ) = B(p 2 ), which in a full treatment with no trun- 
cation would be the case. However, in an incomplete 
treatment where uncontrolled approximations are made 
great care must be taken, because individual diagrams 
contain quadratic divergences that are no longer guaran- 
teed to cancel. 

The simplest solution has been long known and is 
straightforward [32, 33, 58]. In solving the gluon equa- 
tion we only need to determine A or B. As is known 
from perturbation theory, the problematic quadratic di- 
vergences only occur in the A term, so we just project 
onto 2? M „ of Eqs. (15,16) in such a way that only the 



B term is retained. This is achieved using the so-called 
Brown- Pennington projector [59], 



i 



(17) 



where d is the number of dimensions in which we work. 
Using the A term leaves any solution exposed to unphys- 
ical quadratic divergences, and these are present in one 
set of solutions in the literature [8, 57]. The solutions 
are obtained by the insertion of an additional step in 
the iterative procedure, where the quadratic divergence 
is subtracted from the result of the gluon loop in the 
gluon equation. This is not necessarily a safe thing to 
do since in the non-perturbative region the integration 
results are a priori unknown. This point has been made 
previously [23] and is of fundamental importance. 

In the IR analysis which is common to this solution and 
all subsequent refinements, the gluon determined from B 
differs from the gluon determined from A signaling the 
breaking of transversality of the propagator. This may be 
seen in [23] in the IR analysis, where the A and B terms 
are analysed separately and the difference is clear. It may 
also be seen by solving self-consistently for the singular 
solution defined by the A function, and then seeing how 
the solutions change when we switch from the A term and 
the B term. This we show in Fig. 8 and we see a clear dif- 
ference between the two gluon dressings, signalling there 
is an issue with transversality. Moreover the B solution 
is not self-consistent. If the usual iterative procedure is 
followed then no self-consistent solutions can be found 
where the B term satisfies the set of equations [60] . 

We thus determine our gluon from the p^Pv term of the 
propagator which is known to yield the correct logarithms 
in the UV limit. Reference [11] addresses this issue by 
introducing parameters into their vertices to make A and 
B come together. 



A. Self-consistent solutions of both propagators in 
a one-loop only system 

We proceed in the traditional manner by starting with 
the simplest conceivable system, then attempt to obtain 
solutions. In this system, we neglect the two two-loop 
dressed contributions in Fig. 1 to the gluon propagator 
and also drop the quark interaction for now. The re- 
quired input is the triple-gluon vertex and the ingredients 
described above for the ghost equation. 

There are several sources of information regarding the 
triple-gluon vertex. One property that we consider im- 
portant that tends to be neglected is the Bose symmetry 
of the vertex. Bose symmetry is, of course, present at all 
orders in perturbation theory, and in the solution of the 
WSTI for this vertex [6, 36]. 

We proceed by considering three dressings of the 
triple-gluon vertex of increasing complexity, we factor off 
i9f abc , and write the bare vertex Lorentz structure as 
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FIG. 8. (Color Online) The singular ghost solutions and vio- 
lation of transverality in the infrared. The dashed gluon and 
dotted ghost curves are obtained first using the A term and 
then we switch to the B term and we find that transversality 
is broken since the curves differ. The units of p 2 are arbitrary 
internal units. 
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(20) 

The first vertex, T^ A \ has been used to reproduce the 
precise perturbatively resummed one-loop running of the 
gluon propagator [61], note that it does not reproduce 
the one-loop behavior of the triple-gluon vertex itself, 
nor does it have Bose-symmetry between each leg. The 
second vertex, T^ B "> is a simple symmetric vertex inspired 
by WSTI solutions that involve ratios of dressing func- 
tions. Finally, r( c ) is an approximate solution to the 
WSTI itself using a bare ghost-gluon scattering kernel 
(T^ = g^). Since the ghost-gluon vertex WSTI and 
hence the ghost-gluon scattering kernel contributions are 
not precisely known, we do not go further than this at 
present. A full solution of the triple-gluon WSTI using 
an approximate ghost-gluon scattering kernel is available 
in the literature [6, 25, 26], however the associated ghost- 
gluon vertex is insufficient to obtain solutions in a finite 
ghost system. 

We select a conservative set of parameters that from 
experience allow solutions for our vertices given above. 
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These are not intended to describe the physical world 
in any way, but rather to expose the salient features of 
the truncation. For sensible comparison, we only use the 
same parameters, where applicable. 

The solutions obtained are shown in Fig. 9. We label 
the solutions using (i, j) to denote the vertices we use. 

-.(<) 



In order to explain this, we write the renormalised 
(subtracted) gluon equation in the following form, 

GKp 2 )- 1 = GK^r 1 + n 29 b 2 V) + u 2c (p 2 ,^), (21) 

the functions II (p 2 ,/i 2 ) are the loop integrals subtracted 
at the renormalization point /i 2 . These polarization func- 



These are defined above and correspond to for the tions are then zero at p 2 = i± 2 and Q£{y 2 ) = 1. It is 

expected that the gluon propagator dressing should be 
positive for all momenta. 



ghost-gluon vertex and T]il p for the triple-gluon vertex. 
Similarly to the solution of the ghost equation alone, we 
find a strong sensitivity to the vertices we choose, partic- 
uarly in the all-important region that is most relevant to 
physics: roughly between 0.1 and 10 in these momentum- 
squared units (see Figs. (4-12)). Although we do not 
match to a specific scale, the peak of the gluon dressing 
function would be expected to have a close relation to the 
fundamental scale Aqcd- The differences between these 
solutions with the range of vertices {i, j} would undoubt- 
edly have an effect on physical quantities, for example, 
hadron masses and form- factors [62]. 
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FIG. 9. (Color Online) The range of solutions obtained using 
the possible vertex combinations. The label refers to 

the vertices used in obtaining the solutions corresponding to 
Tfj/ for the ghost-gluon vertex and T^} p for the triple-gluon 



vertex. The missing curve corresponds to the and T 
self-consistent solutions were not obtained there. The param- 
eters are not varied between the solutions, only the vertices. 
The units of p 2 are arbitrary since we have not fixed the cou- 
pling to a physical value. 



B. Infrared Loop Contributions from different 
vertices 

We note that previous studies have found ghost- 
dominance in the IR region of the gluon equation. We 
show this to be false and an artifact of choosing a sim- 
ple vertex. We find that as in the perturbative region, 
the ghost and gluon loops contribute a similar numeri- 
cal amount. This may negate arguments relating to a 
ghost-gluon vertex that is transverse alone. 



The ghost-loop polarization function, Il2 C (p , /i ), is 
highly sensitive to the ghost-gluon vertex, and the gluon 
loop polarization function, Il2 ff (p 2 ), is most sensitive to 
the choice of triple-gluon vertex. This is completely as 
expected. There are additional feedback effects that hap- 
pen when the equations are solved together as described 
in appendix A, but these effects are typically smaller. 
These functions diverge as p 2 (or a little less) so we mul- 
tiply up by p 2 and plot the numerical results in Fig. 10. 

In Fig. 10 we see the key feature that we wish to elu- 
cidate, and that is, for the Bose-symmetric triple-gluon 
vertices, the IR gluon-loop contributions are non-zero. 
Both ghost and gluon loops contribute at a similar nu- 
merical order, whilst the ghost must be larger to keep the 
gluon dressing function positive in the IR limit. There 
is clearly some interplay between these, as can be seen 
by the differences between the transverse vertex and the 
modified version that satisfies Taylor's theorem. An in- 
terdependence of this type is to be expected, since the 
ghost-gluon scattering kernel appears separately in the 
WSTI for both vertices, and the idealised full solution 
would uniquely constrain this relation. 

The behavior of these functions is important, partic- 
ularly in light of the above partial cancellation, if the 
gluon propagator dressing function is to be positive ev- 
erywhere. As the loop contribution induced by the triple- 
gluon dressing is negative then some canceling positive 
contribution is required from the loop containing the 
ghost-gluon vertex. This is only present for the two 
vertices we show here. For the bare vertex or a ver- 
tex dressed by simple ratios, no such contribution occurs 
and hence self-consistent solutions are not possible. It 
is this extreme sensitivity, highlighted in Fig. 10, to the 
ansatze for both the ghost-gluon and triple gluon vertices 
that makes us conclude that consistent gluon and ghost 
propagators have not yet been determined in continuum 
strong coupling QCD. 

An alternative to this would be for a positive contri- 
bution to arise from the gluon loop. This would require 
further corrections to be added to the triple-gluon ver- 
tex. Some evidence for this may already exist in recent 
lattice studies [63]. These contributions could be due to 
the WSTI-unconstrained parts of the triple-gluon vertex 
and hence are not present in T^) p of Eq. (20), or a more 
complete WSTI solution [6]. 

A point of note is that the coefficient of p 2 in the IR 
limit of the gluon propagator dressing is determined by 
the sum of the ghost and gluon loop functions in that 
limit. In order to determine this with any precision it is 
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FIG. 10. (Color Online) Gluon polarization functions, p Il2c(p ) (positive values) and p IU^p ) (negative values). Left: IR 
region, Right: UV region. In the UV the curves are indistinguishable on this scale. In the IR, the ghost-loop curves for multiple 
solutions lie on top of each other. The bare ghost-gluon vertex, not shown here, always gives a vanishing IR contribution given 
an IR vanishing gluon dressing and an IR finite ghost dressing. The functions p 2 il2 C (p 2 , fJ 2 ) are labelled according to the input 



vertices with the contributions from r^ 1 ' and T^ p labelled as in the plot. The units of p 2 are arbitrary. 



important to get the vertices right. This coefficient is of 
some importance since it determines the effective gluon 
'mass' term. The mass term, though almost certainly 
gauge dependent, will inevitably affect physical quanti- 
ties like hadron masses and form-factors. Many practical 
models contain mass terms for the gluon; for a recent 
example see [64]. 

Considering these two contributions from the differ- 
ent loops we also find the reason why the gluon dressing 
looks similar in both the infinite and finite ghost solu- 
tions. In the singular solution the bare vertex is used 
and the large contribution from the ghost loop is pro- 
vided by the singular propagator dressing. Conversely 
in the finite ghost solution a non-trivial vertex provides 
a similar large contribution and a qualitatively similar 
vanishing gluon arises. 



V. COMPARISON TO LATTICE QCD 

Primarily we have been motivated by theoretical issues 
encountered in solving the Schwinger-Dyson equations 
for the gluon and ghost propagators. However a com- 
plementary technique is available where these quantities 
have been calculated and that is Lattice QCD. The pos- 
sible issues there are quite different to those that we may 
induce here by truncation, so a comparison is a useful 
independent cross-check. Importantly the lattice com- 
putations, for which there are extremely precise results, 
are for the pure gauge sector we investigate here in the 
continuum. 

Many recent lattice studies exist [65-67] , starting with 
the early work of [2]. We compare our calculations to 
[45] which provides results in Landau gauge for both 
dressing functions. The qualitative behavior there is as 
we have here with a finite ghost dressing function and 
a finite gluon propagator, corresponding to a vanishing 
gluon propagator dressing function. 



In obtaining these solutions we use the preferred set of 
vertices, from Eq. (20) and the ghost-gluon vertex 
given in Eq. (14). We then tune the parameters to ob- 
tain a reasonable representation of the lattice data and 
these are given in Table I. We note several numerical dis- 
crepancies. Most importantly in the UV, the available 
parameter space does not include the coupling strength 
found on the lattice; smaller couplings have to be used. 
The iterative procedure described in Appendix A breaks 
down at larger couplings signalling an absence of solu- 
tions in this given truncation. This is visible in the plots 
as a smaller gradient of the perturbatively resummcd log 
for the SDE solutions. Secondly, the peak of the gluon 
dressing function in Fig. 11 is not as large as found on 
the lattice. It is expected that the two-loop dressings 
would make a contribution here, so that could be the 
source of the difference. Other effects will of course in- 
clude the different coupling values and uncertainties in 
the vertices. We do not comment upon any differences 
due to the method used to extract the predictions from 
the lattice. 



Solution 


A 4 




Ajr 


Mr 


Dotted 


650 


0.1313 


6.0 


3.5 


Dot-dashed 


650 


0.1200 


3.8 


3.7 


Dashed 


650 


0.1225 


7.5 


2.0 



TABLE I. Parameters used in obtaining the tuned self- 
consistent lattice solutions, /i and Air, can be regarded as 
being in GeV units. 



VI. OUTLOOK 

The gluon and ghost propagators are the basic Green's 
functions that embody not only the short distance be- 
havior determined by the asymptotically free nature of 
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FIG. 11. (Color Online) Self-consistent solutions tuned to lattice solutions, showing the dressing functions. The solid curves 
depict a smooth fit to the lattice data. The heavier region is where the functions are represented by lattice data and the 
feint region represents the natural extrapolation. The broken curves are the tuned solutions. Left: The blue peaked curves 
correspond to the gluon dressing function Q£(p 2 ), this vanishes as p 2 — > 0, the red monotonic curves correspond to the ghost 
dressing function Qh(p 2 ). Right: The upper (blue) curves correspond to Q£(p 2 ) /p 2 , this is ~ 10 as p 2 — > 0, the lower (red) 
curves correspond to the ghost dressing function. 



QCD, but confinement dynamics at larger distances. 
This behavior is encoded in solutions of the appropriate 
Schwinger-Dyson equations. Here we first investigated 
the existence of solutions for a simplified ghost equation 
using a fixed gluon input. Multiple solutions were found. 
However applying a condition commonly used in pertur- 
bative analyses, that Qh(n 2 ) — 1 where fi 2 is the point 
in the perturbative region at which both gluon and ghost 
are renormalised, we found just one solution was pre- 
ferred. Using this result we then investigated the neces- 
sary terms required in order to construct a self-consistent 
solution for the coupled gluon and ghost dressing func- 
tions with input interaction vertices. A range of solutions 
resulted that are qualitatively similar to those computed 
on the lattice. However, perfect quantitative agreement 
has yet to be established. This is assumed to be due to 
the approximations made here or on the lattice. We drew 
special attention to differences induced by vertex choices 
and we have argued that the neglected two-loop dressings 
are likely to give rise to significant changes. The hope is 
that constraints on the vertices, particularly that of the 
full ghost-gluon interaction, can be found that uniquely 
specify their structure. At present the ad hoc vertices 
used are motivated more by the practicalities of finding 
solutions than by constraints derived from the fundamen- 
tal field theory. Thus we conclude that consistent gluon 
and ghost propagators have yet to be determined in con- 
tinuum strong coupling QCD. 

In the pure gauge sector studied here, one can define a 
non-perturbative running coupling following Taylor [37] 
from the ghost-gluon vertex renormalization, 

«^ 2 ) = ?^(f 2 )^ 2 ) ( 22 ) 
where /i 2 is the renormalization point and g(n 2 ) fixes the 



physical scale. Experiment with 5 flavors of fermion gives 
a(M|) = 0.118 in the modified minimal subtraction 
(M S) scheme. This can be related to momentum sub- 
traction with zero flavors appropriate to the calculations 
we performed here, by using perturbative results from 
the two schemes [35, 68, 69]. Thus, this value of a in the 
MS scheme with 5 flavors translates into a(M§) = 0.08 
with zero flavors in the momentum subtraction scheme. 
Then with Air = 2 GeV and A/"ir = 2 we obtain the 
function axip 2 ) shown in Fig. 12. 

The behavior of the gluon propagator dressings in 
Figs. 4, 9, 11 and 12 can be interpreted in terms of the 
dynamical generation of an effective gluon mass. This in 
turn leads to the following proposal for a corresponding 
definition of the running coupling 

/ 2n g 2 (M 2 ) P 2 + m 2 (p 2 ) 

CtECiP = — j 2 * h (P )G*\P ' ( 23 ) 

recently developed in a study of the effective charges 
given by dressing functions with such qualitative be- 
havior [55]. Applying this to the solutions found using 
a(M 2 ) = 0.08, Air, = 2 GeV and Mr = 2 we find 
m g (0) = 350 MeV and the function a E c{p 2 ) in Fig. 12. 
The result is not particularly sensitive to the precise form 
of the function m g {p 2 ), however we use that given in [55] 
which is suppressed in the UV. 

Both gluon and ghost dressings do not rise as steeply as 
the lattice-like solutions of Fig. 11. Nevertheless, the col- 
lective effect of the gluon dressing and this running cou- 
pling would just have sufficient strength to induce the dy- 
namical generation of quark mass, cf. the model of Maris 
and Tandy [62, 70]. The inclusion of quarks is, of course, 
an important aspect of making connection to physical 
quantities, like the pion mass and its decay constant. A 
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FIG. 12. (Color Online) The running coupling ar(p 2 ) from 
Eq. (22) is shown as the black solid curve. The running cou- 
pling otEc{p 2 ) from Eq. (23) is shown as the green solid curve. 
The gluon dressing Q£(p 2 ) is the dashed blue curve and the 
dotted red curve is the ghost dressing function Qh(p 2 ). The 
units of p 2 are now very close to GeV 2 . 



complete study requires a simultaneous solution of the 
quark SDE with those of the gluon and ghost. Research 
over the past 15 years indicates that we know what the 
quark propagator functions look like for a whole range 
of quark masses (and hence flavors) with model gluon 
inputs. Indeed modeling such ghost and gluon propaga- 
tors and their interactions, and inserting these into the 
Bethe-Salpeter equations to compute hadronic observ- 
ables is now well advanced [62, 64, 70, 71]. However, a 
comprehensive computation with a self-consistent inves- 
tigation of the coupled quark, gluon and ghost equations 
in continuum QCD, followed by the development of a 
trustworthy hadron phenomenology, is still in the future. 
The issues exposed by the present study show we have 
some way to go before we can claim robust results from 
such an ab initio approach. 
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Appendix A: Numerical Method 

The equations are solved in Euclidean space using the 
well known method set out in [72]. Using the standard 
spherical four-dimensional coordinate system the two in- 
ner angles are symmetric for the propagator integrations 
so may be integrated out. We are then left with the out- 
ermost angle which we define as the angle between the 
incoming propagator momentum p and the loop coordi- 
nate £. 

The integrals are logarithmically divergent and are reg- 
ularised by integrating to some momentum cutoff n. This 
appears in the renormalization procedure. However the 
results do not depend on it. The functions are subtracted 
at /i 2 and this point is also used to match on to the per- 
turbative result. That is, above this point the one-loop- 
resummed perturbative result is used as an extrapolation 
when required by the loop integrations. 

The functions are represented by Chebychev polyno- 
mials. These are particularly useful for both their in- 
terpolation properties and for the iterative procedure. 
They are mapped onto a logarithmic scale typically in 
the range p 2 g [lO -8 , 10 5 ] although the precise numbers 
are not important. The last Chebychev zero is mapped 
to the point p? since this point is always exactly repro- 
duced by the interpolation. In the IR then the functions 
have to be represented by some extrapolation. In finding 
these finite solutions a constant is sufficient for the ghost 
dressing function and the c g p 2 term is appropriate for the 
gluon dressing, where c g is calculated to match smoothly 
onto the lowest point represented by the Chebychev poly- 
nomials. 

In performing the iterations we use both the Newton- 
Raphson and natural iterative procedures. For the finite 
ghost solutions, then natural iterative procedures are ad- 
equate for finding solutions in both equations. A good 
starting point is always useful and often a requirement, 
particularly at large couplings. We use Eq. (13) for the 
gluon starting point and Qh(p 2 ) = 1 is a sufficient start- 
ing point for the finite ghost dressing. 



Appendix B: Formulae 
1. The Ghost Equation 

The renormalised ghost equation used is, 

GHp 2 )- 1 = Z 3 (^ 2 , K 2 ) +g 2 N c J 

a(-£_, P ,i + ) ic( P ,e+,£-) gh(£+)ge(t-) (bi) 

where £± = £ ±p/2 where £ is the loop integration mo- 
mentum, Z-j, is the ghost renormalization constant. The 
function a multiplies the term in the ghost-gluon ver- 
tex, as defined in Eq. (6). The function K, arises from the 



tensor contractions in the Feynman rules and is given by, 



moves the term Z 3 . 
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£ 2 sin 2 9 

where 9 is the integration angle defined via £.p — 
\£\\p\cos9. This equation is typically subtracted from 
itself at p? which is the renormalization point. This re- 



2. The Gluon Equation 

As described in the text the gluon equation is made 
up of more terms and contains two dressed one-loop in- 
tegrations in the formulation we have used here. We give 
no details of the tadpole term since it yields results that 
are proportional to the g^ v term of the propagator and 
hence does not contribute with the method we describe 
above. The gluon equation that we have used is given by 



giiP 2 )- 1 =z 3 (h\k 2 ) + n 2 > 2 , k 2 ) + n 2c ( P 2 , k 2 ) (B3) 
=QI(p 2 Y 1 + n 2g ( P 2 , k 2 ) - n 2g (p 2 lK 2 ) + n 2c ( P 2 , « 2 ) - n 2c ( M 2 , K 2 ) (B4) 
=Qi(u?y x + n s 2 f (p 2 , fj, 2 ) + n^V, m 2 ). (B5) 



where in the second line we subtract the equation from 
itself at the point p 2 and in the third line we introduce the 
subtracted loop integrations IFj ub (p 2 , p 2 ) = Hj(p 2 7 n 2 ) — 
Tlj(n 2 , k 2 ) which for each diagram are actually all the 
same in a properly renormalised system when everything 
above the point p 2 is assumed to be described by one-loop 
perturbation theory so we drop the (sub) notation below. 
We next define the content of the loop integrations in the 



II functions. We give this only for the most complicated 
case: of vertices and r( c ) since the others may be 
straightforwardly deduced from these by setting factors 
to 1 and/or multiplying by the appropriate dressings of 
the vertex. 

First we give the ghost-loop contribution. This con- 
tains the two functions a(k,p,q) and /3(k,p,q) from the 
ghost-gluon vertex dressings of Eq. (6) that differ be- 
tween and 



n 2c ( P 2 V)=(-i) 
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where d — 4 is the number of dimensions, g is the value of the coupling at the renormalization point, the functions a 
and (3 arise from the two terms in the ghost-gluon vertex, Eq. (6). The kinematic terms are contracted together into 
the two M functions, which are derived to be, 



M a (p,£) = 
M p (p,£) = (l-d)p.£_ 
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pr_ 

4 



(B7) 
(B8) 



For the bare vertex, the (3 term is zero and for the transverse vertex, T^, then a(k,p, q) — 1. For the Taylor-respecting 
vertex r^ 2 ' then we replace p^ — > fc p + q^ in Eq. (14) and read off the respective coefficients. Similarly for the gluon 
loop we have, 



n 29 b 2 ) = 



N c g 2 ( P 2 
(d-1) 



d 4 £ Q£(£\)Q£{£ 2 _ 



where the ratios of dressing functions are chosen to be 

oi\ Qh(£ 2 _^ 1 /r.hfrfi 



_ _ l f Gh(f) 



G£(p 2 



112 2 \QW + 



(2tt) 4 4(pH 2 + £ 2 _) 2 
GKp 2 ) 



(niQi + n 2 Q 2 + n 3 Q 3 ^ (B9) 



G£{£ 2 



_1 (Qh(l\) Qh{l\), 

n "-2\gmr) + w) ] (B10) 
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and the associated momentum contractions are, 

Qi - 2{£.p 2 - £ 2 p 2 ) [4£ 2 (5p 2 + 2d£.p) + p 2 {7p 2 + 2(12 - d)Lpj\ , (Bll) 

Q 2 = -Z2d{£.p) 4 + (3d - 4)/(£.p) 2 - 4S£ 6 p 2 + 24£ 4 (2d£.p 2 - p 4 ) + £ 2 p 2 (8(d + 6)£.p 2 + p 4 ) , (B12) 

Q 3 = 2(£p 2 - £V) [4^ 2 (5p 2 - 2d£.p) + p 2 (7 P 2 - 2(12 - d)Lp)] . (B13) 



In order to obtain a vertex proportional to the bare 
Lorentz structure as in r( A ) and r( B ) then we just have 
to make a replacement, for example for we have, 



Q£(£\)Q£(£ 2 _) ' 



(B14) 



for all three TZi's, and a similar substitution is made to 
obtain the equation for the dressing. 

The smoothed step function used for interpolating be- 
tween IR and UV vertex functions is taken from [11]. 
However, a range of forms may be used, 



A IR 



(fc 2 + A 2 R )(^ + Al 2 R )( g 2 + A 2 R ) 



(B15) 



where Air gives the change-over between the perturba- 
tive and non-perturbative form. The solutions are sen- 
sitive to the choice of this parameter. In the modified 
form of the vertex of Eq. (14) an extra parameter A/ir 
multiplying the non-perturbative term was introduced. 



3. Resummed One-loop running 

In the UV we match the gluon and ghost propagator 
dressings on to the resummed one-loop form from per- 
turbation theory. This is necessary as some extrapola- 
tion beyond the Chebychev region in the UV is always 



required. The form that we find to be most useful nu- 
merically is given by, 



gi(p 2 ) = g£(n 2 )(a^ 2 )/a( P 2 )) d ° 

gh(p 2 )=gh(p 2 )(a(n 2 )/a(p 2 )) d c 



(B16) 
(B17) 



where d g = —13/22 (cf. the modeling in Eq. (13)) and 
d c = —9/44 (as appropriate for the pure gauge sector) 
which is applied where p 2 > n 2 . Numerically, we also use 
the standard relation for the running coupling, 



1 



1 



a(p 2 ) ct(^ 2 ) 



i UN C p 2 

127T [l z 



(B18) 



where ct(ji 2 ) is specified to fix the scale. 

Since only the triple-gluon vertex Tffip perfectly repro- 
duces the perturbativcly resummed one-loop running, it 
is useful to interpolate between this in UV and the other 
vertices in the mid-momentum region and below. The 

(C) 

WSTI triple-gluon vertex, Tp Vp in particular induces a 
steeper gluon dressing function than would be expected 
from perturbative studies. Some deviation is allowed due 
to the effects of higher orders pcrturbatively and since we 
are still at the modeling stage then interpolating between 
these two vertices between an IR scale A 2 R and ^i 2 ~ M§ 
is currently the best that can be done. The transition is 
implemented here using Jir (see Eq. (10) and Eq. (B15)) 
and (1 — Jir) to multiply the respective terms. 
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